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Introduction 


One of the primary goals of NASA’s high-speed research program is to 
determine the feasibility of designing an environmentally safe commercial su- 
personic transport airplane. The largest environmental concern is focused 
on the amount of ozone destroying nitrogen oxides (NO*.) that would be 
injected into the lower stratosphere during the cruise portion of the flight. 
The limitations placed on NO x emission require more than an order of mag- 
nitude reduction over current engine designs. To develop strategies to meet 
this goal requires first gaining a fundamental understanding of the combus- 
tion chemistry. 

To accurately model the combustor requires a computational fluid dy- 
namics approach that includes both turbulence and chemistry. Since many 
of the important chemical processes in this regime involve highly reactive 
radicals, an experimental determination of the required thermodynamic data 
and rate constants is often very difficult. Unlike experimental approaches, 
theoretical methods are as applicable to highly reactive species as stable 
ones. Also our approximation of treating the dynamics classically becomes 
more accurate with increasing temperature. In this article we review recent 
progress in generating thermodynamic properties and rate constants that are 
required to understand NO x formation in the combustion process. We also 
describe our one-dimensional modeling efforts to validate an NH3 combus- 
tion reaction mechanism. We have been working in collaboration with Marty 
Rabinowitz at Lewis research center, to ensure that our theoretical work is 
focused on the most important thermodynamic quantities and rate constants 
required in the chemical data base. 
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Outline 


Thermodynamic properties 

1. C-H bond dissociation energies of CH4 

2 . Resolution of the controversy concerning the C-H bond energy of 
acetylene 

Reaction rate constants 

1 . Lifetime of HN 2 to unimolecular decay 

2 . C 2 H 2 + C 2 H 2 -> ? 

3. propane + O — ► propyl + OH 

Modeling 

1. ID modeling to validate an NH 3 combustion mechanism (H-N-O) 

2 . ID modeling of turbine combustors with CH 4 fuel to identify impor- 
tant H-C-N-O reactions. 



Calculation of heats of formation 


Although the heat of formation is well known for many molecules, indi- 
vidual bond energies are often very uncertain. We have developed an inex- 
pensive computational approach for predicting accurate (1 to 2 kcal/mole) 
C-C and C-H bond energies [1-2]. Currently we are attempting to assess the 
accuracy with which we can compute N-H, N-N, O-H and 0-0 bond ener- 
gies [3]. An accurate knowledge of bond energies is required to model the 
chemistry occuring in jet engines. For purposes of calibration, we have stud- 
ied the successive bond energies of methane [1], These are compared with 
experiment in the following table. At a modest level of theory, our directly 
calculated values are about 2 kcal/mole too small after the inclusion of the 
vibrational zero-point energy (ZPT). However, 0.5 kcal/mole accuracy is 
achievable with very large calculations. When we add 2 kcal/mole to correct 
for deficiencies in the theoretical bond energies (column labeled estimate), 
excellent agreement is obtained with experiment. This gives us consider- 
able confidence that our results for hydrocarbons and other systems, such as 
ketene [2] and methanol [4] are also very accurate. 
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C-H dissociation energy in C 2 H 2 


There is currently considerable controversy regarding the C-H disso- 
ciation energy of acetylene, DfJ(HCC-H). This is in large part due to the 
low and presumably very accurate value determined by Green, Kinsey 
and Field (GKF) [5] using Stark anticrossing spectroscopy. Their result, 
126.647(2) kcal/mole, is claimed to be an upper bound , and is consistent with 
the 127±1.5 kcal/mole estimate of Segall et al. [6] obtained by measuring the 
kinetic energy (K.E.) of the hydrogen atom fragment using Doppler mul- 
tiphoton ionization spectroscopy. However, these values are substantially 
smaller than previous theoretical estimates [7,8] and other recent experimen- 
tal results, such as the Do value of 131.3±0.7 kcal/mole measured by Ervin 
et al. [9] using the techniques of negative ion photoelectron spectroscopy and 
gas phase proton transfer kinetics. 

We undertook a systematic study of the C-H bond dissociation energy 
in acetylene with respect to improvements in the level of theoretical treat- 
ment [10]. Our best estimate for the D 0 value of 130.1±1.0 kcal/mole is in 
good agreement with the recent experiment of Ervin et al. [9] and previous 
theoretical calculations [7,8]. We believe our error bars to have a better than 
90% probability of bracketing the correct D 0 value, and we therefore seri- 
ously question the recent upper bound inferred from Stark anticrossing spec- 
troscopy [5]. Detailed theoretical studies of the C 2 H 2 and C 2 H vibrational 
frequencies have also led us to revise upward the D® values determined from 
kinetic data. 


<126.647 ±0.002 
127. ±1.5 
124-127 
13 1.3±0.7 
132. ±2 
132.6±1.2 


130.1±1.0 


Experimental 

Stark anti-crossing spectroscopy 
HCCH + hu -» C 2 H + H(K.E.) 
Kinetic data 

AH QC!d (HCCH) + EA(C 2 H) - IP(H) 
HCCH + hu -+ C 2 H(K.E.) + H 
photoionization of HCCH 

Theoretical 

MRCI 
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Reaction Rate Constants 


The first step in determining an accurate rate constant is to compute 
the potential energy surface — that is, for “all” possible arrangements of the 
atoms, solve for the electronic energy. The surface must then be fit to an 
analytic form in order to solve for the nuclear motion on this surface using 
either classical or quantum mechanical scattering methods We are currently 
exploring the use of variational transition state approaches with tunneling 
corrections for determining qualitative (factor of 2-3) reaction rate constants 
and product branching ratios. One of the advantages of this method is that 
it requires only a small portion of the complete potential energy surface. 
This approach has been applied to determining the rate constant of the re- 
action of propane with atomic oxygen — see later discussion. 

Several reactions are currently under study in our laboratory. For exam- 
ple, the reaction 

H + 0 2 -> HO + 0 

is being studied, since this reaction consumes most of the 0 2 in typical hy- 
drocarbon combustion. However, the rate constant for this reaction is un- 
certain by a factor of six at flame temperatures. A global potential energy 
surface(PES) has been developed and fitting is in progress [11-14]. 

Another molecule that we have studied extensively is HN 2 , since it has 
been postulated as an important species in thermal De-NO*. processes. A 
global PES has been developed for 

H + N 2 -> HN 2 

and an accurate lifetime has been computed for the HN 2 molecule [15-18]. 

The reaction of C 2 H 2 with itself is postulated to be important in the fi- 
nal steps of hydrocarbon combustion. We believe that acetylene must first 
be converted to the vinylidene isomer (CCH 2 ) before reaction can occur — 
see later discussion. 

The reaction 

GH + N 2 -> HCN + N 

is believed to be the rate determining step in “prompt NO z ” formation. A 
large number of stationary points have been located on this potential energy 
surface. Large-scale investigations of the surfaces are in progress. Finally, 
work is beginning on the CH 3 +0 2 and CH 3 +OH reactions to determine 
product branching ratios. Model studies have shown these to be important 
in the combustion of jet fuel. 
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Global potential energy surface for H+N 2 — + HN 2 

The HN 2 species has been postulated as an important intermediate in 
the thermal De-NO r process. The HN 2 species was found to be unstable 
with respect to H+N 2 by 3.0 kcal/mol, but quasibound by 12.2 kcal/mol 
due to a barrier to dissociation — see figure below. The computed N 2 -H life- 
time [17] (based on 1-D tunneling through an Eckart barrier) was five orders 
of magnitude smaller than the value assumed in the experimental analysis. 
Koizumi and Schatz (Northwestern) have fit this surface and carried out 
coupled channel calculations to determine the lifetime. These rigorous cal- 
culations confirm those based on a Wigner model of tunneling that indicate 
this radical has a very short lifetime to unimolecular decomposition [18]. 

This suggests that commonly used reaction mechanisms for NO* chem- 
istry such as the Sandia reaction set [19] are incorrect. 
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c 2 h 2 +c 2 h 2 ? 

The reaction of C 2 H 2 with itself is postulated to be important in the fi- 
nal steps of hydrocarbon combustion. It is exceedingly unlikely that any re- 
action takes place between two acetylene molecules at the temperatures of 
interest, but the vinylidene (CCH 2 ) isomer is only about 30 kcal/mol above 
acetylene, and a number of products can result from reaction of vinylidene 
with acetylene. 

So far, seven stable isomers of C 4 H 4 have been found using self- 
consistent-field (SCF) and multiconfigurational self-consistent-field (MC- 
SCF) wave functions, as well as several stationary points that are not min- 
ima or are not stable with respect to symmetry-breaking nuclear distortions. 
The lowest energy C 4 H 4 isomer is vinylacetylene, which is about 45 kcal/mol 
more stable than two acetylenes. Butatriene is about 10 kcal/mol higher in 
energy, and methylenecyclopropene is about 14 kcal/mol higher still. An- 
other isomer, cyclobutadiene, is only a little more stable than two acetylenes. 
The other stationary points are much higher in energy. The energetics here 
are not fully refined at the correlated level, but correlation does not seem to 
have a significant effect on the thermochemistry. 

Several pathways have been followed on this surface using large MCSCF 
configurational spaces, in preliminary attempts to find transition states for 
the first reaction steps of vinylidene with acetylene, but while these have 
yielded another fairly stable isomer, carbenecyclopropane, and transition 
states for some rearrangements have been found, the initial transition states 
have not yet been located. However, the investigation is at an early stage 
and there should not be a problem determining the reaction pathways. 


974 



CM 3 MOLECULE" 


The Reaction of Propane with Atomic Oxygen 

We have studied the hydrogen abstraction from the central carbon of 
propane by ground state oxygen [20]. The energy barrier on the potential 
energy surface was determined using SCF gradient methods. A more accu- 
rate barrier and exoergicity were determined using correlated methods. Fur- 
thermore, a correction factor to account for deficiencies in the basis set and 
electron correlation treatment was determined by carrying out large scale 
calculations for the analogous reaction H 2 + O — ► H -f- OH. This correction 
lowers the saddle point and product energies by 4.00 kcal/mol relative to the 
reactants. The final computed energies are 4.30 kcal/mol for the barrier and 
7.00 for the exoergicity. These compare favorably with the experimentally 
determined values of 5.0 and 7.0 kcal/mol. 

The potential energy profile, molecular geometries and vibrational fre- 
quencies were used for a transition state theory calculation of the rate con- 
stant for the hydrogen abstraction reaction. Two approximate treatments 
of tunneling were carried out. The Wigner model assumes that the energy 
barrier is an inverted parabola truncated at the reactant and product en- 
ergy limits. The Eckart model uses a potential energy function based on a 
modified hyperbolic secant. It can be seen from the figure that the transi- 
tion state theory calculations with either tunneling model agree with the ex- 
perimental data to within a factor of two over the entire temperature range 
(300-2000K). Greater than a factor of two accuracy in the rate would require 
a global surface, which is exceedingly difficult for a system with this many 
atoms. 





Thermal de-NO x 

Experiment of Lyon and Hardy [21] showing NO removal by ammonia 
injection. The process works in a narrow temperature range, and oxygen 
must be present for NH 2 formation. Initial conditions of experiment are: 
NO=225 ppm, NH3=450 ppm, 0 2 =1.23%, remainder He. The pressure is 
1.1 atm and the observations are made at 0.1 sec. Case two is similar to 
case one except for 225 ppm H 2 displacing an equal amount of He. Modeling 
shown here uses the reaction mechanism of Miller and Bowman [19], with a 
long lifetime for N 2 H (r=10 -4 sec). Calculations assume a plug flow reactor, 
and are simialr to Miller and Bowman [19]. The same calculations done with 
a shorter N 2 H lifetime, in an attempt to match the theoretical value of the 
lifetime (< 10~ 9 sec), are unable to reproduce the experimental results. 
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Conclusions 


Computational chemistry techniques have utility in computing accu- 
rate thermodynamic properties such as C-H bond energies. Accurate C-H 
bond energies have been computed for methane, ethylene, ketene, acetylene, 
and methanol. Reaction rates can be computed at least to within a factor of 
two. The reactions currently under study in our laboratory include H+ 0 2 — > 
HO-fO, H+N 2 -> HN 2 , C 2 H 2 with itself, CH+N 2 HCN+N, CH 3 +OH, and 
CH3+0 2 . Our calculations demonstrate that the HN 2 molecule has a short 
lifetime (<10 -9 sec) to unimolecular decomposition. We are presently trying 
to incorporate this fact into model studies of the effect of NH3 on NO re- 
moval. We have been working in close collaboration with Marty Rabinowitz 
at Lewis to ensure that we are addressing the key chemical issues involved in 
reducing NO emission in jet engines. 
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